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1.  INTRODUCTION 


In  the  area  of  signal  processing,  eigenstructure  methods  are  becoming  popular  in 
determination  of  the  number  of  signals  and  estimation  of  the  parameters  of  the 
when  noise  is  present  Problems  of  signal  detection  in  presence  of  noise  are 
complex  and  they  arise  in  many  situations.  For  example,  in  the  area  of  active 
detection,  signals  are  transmitted  and  the  echoes  from  the  targets  are  recorded 
observations  are  corrupted  by  disturbances  in  the  sea  and  we  are  interested  in 
identification  of  the  targets.  In  radar  detection,  pulses  are  sent  and  the  reflections  from 
the  targets  are  recorded  Here,  the  observations  may  be  subject  to  atmospheric 

disturbance  and  disturbances  due  to  jammer  and  clutter.  Eigenstructure  methods  were  used 
in  the  past  by  Haykins  (1985),  Kumaresan  and  Tufts  (1980),  Ligget  (1973),  Owsley 
((1977),(1985)).  Pisarenko  (1973),  Redd!  (1979),  Schmidt  (1979),  Tufts  and  Kumaresan  (1982) 
and  Wax,  Shan  and  Kailath  (1984).  Recently,  Wax  and  Kailath  ((1984),(1985)).  Zhao, 
Krishnaiah  and  Bai  ((1985a),  (1985b)),  Bai,  Krishnaiah  and  Zliao  (1985)  and  Paulraj  and  Kailath 
(1986)  considered  the  problem  of  determination  of  the  number  of  signals  using  nradel 
selection  approach.  The  object  of  this  paper  is  to  give  a  review  of  the  work  on  model 
selection  approach  to  determine  the  number  of  signals  in  presence  of  noise. 

In  the  model  considered  in  this  paper,  the  observations  from  the  sensors  are 
expressed  as  the  sum  of  noise  and  linear  combinations  of  the  wavefronts  of  the  signals 
from  the  sources^  When  the  noise  is  spatially  white,  the  number  of  signals  is  related  to 
the  multiplicity  of  the  smallest  eigenvalue  of  the  covariance  matrix  of  the  observation 
vectors.  Various  test  procedures  are  available  in  the  literature  (e.g,  see  Anderson  (1963). 
Bartlett  (1954),  Krishnaiah  (1976)  and  Rao  (1973))  for  testing  the  hypothesis  that  the 
multiplicity  of  the  smallest  eigenvalue  of  the  covariance  nutrix  is  equal  to  a  specified  value 
when  the  underlying  distribution  is  multivariate  normal  Tfie  above  tests  can  be  used  to  test 
the  hypothesis  that  the  number  of  signals  is  equal  to  a  specified  value  Wfwn  the  noise  is 
spatially  colored,  the  number  of  signals  is  related  to  the  multiplicity^  of  the  smallest 
eigenvalue  of  Z.Z  ^where  1  is  the  covariance  matrix  of  the  noise  vector.  Rao  (1983) 
considered  the  following  related  problem  Let  and  be  distributed  independently  as 
central  Wishart  matrices  with  n  and  n  degrees  of  freedom  and  E(S)  =  nZ  for  i  =  1,2. 

12  i  i  i 

Also,  lot  Z„  =  r  XZ,  where  X  is  a  known  or  unknown  scalar,  and  T  is  nonnegative 


sisi 


definite  metrix.  Rao  (1983)  derived  the  likelihood  ratio  test  statistic  for  testing  the 
hypothesis  on  the  rank  of  r  when  X  is  unknown;  when  X  is  known,  he  derived  modified 
likelihood  ratio  test  statistic  The  results  of  Rao  can  be  used  to  test  the  hypothesis  that 
the  number  of  signals  is  equal  to  a  specified  value  in  the  case  of  spatially  colored  noise 


There  are  a  number  of  situations  whore  the  experimenter  is  confronted  with  the 
problem  of  making  one  of  the  various  possible  decisions  instead  of  testing  one  specified 
hypothesie  ki  these  situations,  one  can  formulate  the  problem  within  the  framework  of 
model  selectioa  In  our  case,  we  have  to  choose  one  of  the  several  possible  models 
where  each  model  corresponds  to  the  statement  that  the  number  of  signals  is  equal  to  a 
particular  value.  In  the  model  selection  approaches  discussed  in  this  paper,  we  do  not 
need  any  threshold  value  (critical  value).  In  tests  of  hypotheses,  we  need  a  critical  value 
for  implementation  of  the  tests  and  it  is  quite  complicated  to  compute  the  critical  values  in 
some  situations.  The  model  selection  approaches  discussed  here  involve  using  certain 
information  theoretic  criteria.  Of  course,  model  selection  methods  should  have  optimum 
properties  to  have  an  appeal  One  such  desirable  property  in  our  case  is  the  consistency 
of  the  estimate  of  the  number  of  signals.  Consistency  property  of  the  model  selection 
methods  is  also  discussed  in  this  paper.  A  description  of  the  organization  of  the  paper  is 
given  below. 


In  Section  2,  we  give  some  preliminaries.  In  Section  3,  we  give  a  review  of  the 
literature  on  detection  of  the  number  of  signals  in  presence  of  white  noise.  Analogous 
problems  are  reviewed  in  Section  4  when  the  noise  is  spatially  colored.  In  Section  5,  we 
review  the  literature  on  detection  of  the  number  of  signals  when  the  subarrays  are  located 
at  wide  distances  apart  The  main  emphasis  of  this  review  is  on  the  Efficient  Detection 
Criteria  (EDO  procedures  developed  by  Zhao.  Krishnaiah  and  Bai  ((1985a),  (l98Sb))  very 
recently.  These  procedures  are  strongly  consistent  in  the  case  of  white  noise,  Bai 
Krishnaiah  and  Zhao  (19  SjT a)  gave  exponential  type  bounds  on  the  probability  of  wrong 
detection  of  the  number  of  signals  if  the  EOC  procedures  are  used  These  bounds  tend  to 
zero  rapidly  as  the  sample  size  increases.  Sknilar  bounds  can  be  derived  in  the  case  of 
spatialiy  colored  noise. 


/  IT J  IT-  rj  itrj  rf--  *r^  try  w'JV^  W  <r.^/H 
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2.  PRB.IM1NARIES 

We  now  define  complex  multivariate  normal  and  Wishart  distributions  since  these 
definitions  are  needed  in  the  sequel 

Let  X  =  y  -«■  iy  where  (y*  V  )  is  distributed  as  2p-variate  normal  with  mean  vector 
-  i2  i2 

0  and  covariance  matrix 


j  =  /-I  and  and  are  of  order  p  x  p.  Then  x  :  p  x  1  is  known  to  be  distributed 
as  complex  multivariate  normal  with  mean  vector  0  and  covariance  matrix  £  =  2(£^-fi£2). 
The  density  of  x  is  of*  the  form 

fW  =  {l/if’*|£l}exp(-i£“’x»  (2.2) 

where  a  and  A  respectively  denote  the  complex  conjugate  of  a  and  transpose  of  A  Next 
let  X  .„.x  be  distributed  independently  and  identically  as  complex  multivariate  normal  with 

n 

mean  vector  0  and  covariance  matrix  £.  Than  S  «  7  x.<  is  known  to  be  distributed  as 

j., 

the  central  complex  Wishart  distribution  with  n  degrees  of  freedom  and  E(S)  »  n£.  The 
complex  multivariate  normal  distribution  was  considered  by  Wooding  (1956)  whereas  the 
complex  Wishart  distribution  was  investigated  in  Goodman  (1963).  Complex  multivariate 
normal  distribution  belongs  to  the  family  of  complex  elliptically  symmetric  distributions 
studied  by  Krishnaiah  and  Lin  (1984).  For  a  review  of  the  literature  on  complex  multivariate 
distributions,  the  reader  is  referred  to  Krishnaiah  (1976). 


a  DETECTION  OF  NUMBER  OF  SIGNALS  IN  PRESB^  OF  WHITE  NOISE 

In  this  section,  we  discuss  certain  methods  of  detection  of  number  of  signals  in 
presence  of  white  noise.  The  model  that  is  considered  is  given  by 


(3.1) 


x(t)  =  As(t)  +  n(t) 


5 


vwhere  A  =  [A(*U,A(*)],  xW  =  (x  wr.  sW  *  ($  It)),  and  n'lt)  =  In  (t)._.n  It)). 

-1-q-  Ip  -  1q  -  ip 

liere  the  ataments  of  j-t)i  column  of  A  dapand  upon  thia  paramatars  of  j-th  signal  Also, 
mt)  and  sit)  ars  distributad  indapandantly  as  complax  multivariate  normal  with  Enit)  »  Q 

=  Q  E{n|t)nW')  =  o^l .  E{s|t)si5'}  =  ^».  The  covariance  matrix  I  of  is  given  by 

WWW  Dww  2* 


1,  =  Ai^A  ♦  ah 

2  P 


(3.2) 


In  the  model  13.1),  x|t)  denotes  the  observation  on  j-th  sansor,  n.lt)  denotes  the  noise 
i  J 

component  associated  with  j-th  sansor,  and  is  the  signal  due  to  j-th  wavefront  Also, 

we  assume  that  a  uniform  linear  array  of  p  identical  sensors  is  used  and  q  is  unknown  In 

addition,  we  assume  that  x|t  ),_,x|t )  are  n  independent  observations.  One  special  case  le.g, 

•  1  •  n 

see  Shan,  Wax  and  Kailath  11985))  of  the  A  matrix  is 

Al^)  =  (1,exp|-iw  T  ),-,exp|-i«  lp-1)T  ))  (3.3) 

-j  ^  0  j  0  j 

where  T  =  Id  sinO  )/c,  d  is  the  distance  between  the  sensors,  c  is  the  propagation  speed 

)  J 

of  wavefront  and  0^  denotes  the  angle  of  arrivai  of  j-th  signal 

One  important  problem  in  signal  processing  is  to  estimate  q,  the  number  of  unknown 
signale  We  will  now  discuss  this  problem  from  the  point  of  view  of  model  sdectioa 

Let  X  ^  X  denote  the  eigenvalues  ofZandieti.> _ >_  I  denote  the 

Ip  1  p 

•  •  "  _ 

eigenvalues  of  I  where  n£  =  7  x|t)xtt)’.  Here  we  assume  ttiat  xlt  )_,xlt )  are  distributed 

.  i-  I  -  1  .  n 


2  2 

independent  of  each  other.  Now,  let  H  :  X^  >  X,^,  «  -  «  X  =  o  where  o  is 
unknown  Tlien  denotes  the  hypothesis  that  the  number  of  signals  is  equal  to  k.  It  is 
known  that  the  likelihood  ratio  test  statistic  for  testing  is  given  by 


s  A, 


(3.4) 


Since  ~2logL^  is  distributed  asymptotically  as  Chi-square  with  ip-kHp-k-i-D-l  degrees  of 
freedom,  we  can  use  the  above  statistic  for  testing  for  large  samples.  In  most  of  the 
situations,  we  do  not  know  as  for  what  values  of  k  we  should  test  Also,  if  is 
rejected,  we  should  test  for  conditional  upon  the  rejection  of  H^.  The  distribution 


problem  in  this  case  is  complicated  We  now  discuss  a  model  selection  approach 
proposed  by  Zhao.  Krishnaiah  and  Bai  (1985a)  for  determination  of  the  number  of  signals. 


Let 

EDC(k.C  )  =  -logL^  +  C  (3.5) 

n  k  n 

where  ~  M2p-lc-»-1)-i-1  is  the  number  of  free  parameters  that  have  to  be  estimated 

under  H  and  C  is  chosen  such  that 

k  n 

lim  {C  /n>  =  0  (3.6) 

n^*  " 

lim  {C  /logiogn}  =  •.  (3.7) 

-  n 


Also,  let 

EDCI^C  )  »  min (eSaO.C  )....B3C(p- 1  .C  )} .  (3.8) 

fi  n  n 

A 

Zhao.  Krishnaidh  and  Bai  (1985a)  showed  that  q  is  a  strongly  consistent  estimate  of  q.  TNs 
criterion  is  known  to  be  Efficient  Detection  Criterion  (BDQ.  When  C  -  1  in  (3^5). 

n 

A 

mC(4C  )  reduces  to  Akaike's  AlC  Criterion  but  C  does  not  satisfy  condition  (3.7)  in  this 
n  n 

_  A 

case  When  C  =  1/2logn,  EDC{<^C  )  reduces  to  Schwartz-^issanen’s  minimum  description 
n  n 

length  (MOL)  criterion  and  it  satisfies  conditions  (3.6)  and  (3.7).  Wax  and  Kailath  (1984) 
used  AlC  and  MDL  critaria  for  determination  of  the  number  of  signals.  They  have  pointed 
out  that  the  AlC  criterion  is  not  consistent  whereas  the  MDL  criterion  is  consistent;  their 
proofs  are  incomplete  (see  Zhao.  Krishnaiah  and  Bai  11985c)).  Bai,  Krishnaiah  and  Zhao 
(1985)  gave  an  exponential  bound,  under  certain  conditions,  on  the  probability  of  wrong 
detection  and  this  bound  tends  to  zero  rapidly  as  n^*. 

Let  H*  :  X  >  ~  =  X  ==  1.  Zhao.  Krishnaiah  and  Bai  (1985a)  derived  the 

k  k  k-«'1  p 

following  expression  for  the  likelihood  ratio  test  statistic  L*  for  H*; 

,  P 

logL  =  n  7  (log£  +  1  -  ^)  (3.9) 

k  .  i  i 

i»  1  ♦minftk) 


where  t  denotes  the  number  of  eigenvalues  i  which  are  greater  than  one.  They  also 

A  A 

showed  that  q  is  a  strongly  consistent  estimate  of  q  where  q  is  chosen  such  that 

HDC(q,C^)  =  -logL*  +  (3.10) 

where  satisfies  the  conditions  (3.6)  and  (3.7)  and  y^>P>  =  k(2p-k-»-1). 

In  the  S)C  procedure  described  above,  it  was  assumed  that  x(t  ),-,x(t )  are 

•  1  »  n 

distributed  independently  as  complex  multivariate  Gaussiaa  But  these  assumptions  need  not 

A 

be  valid  in  a  number  of  situations.  But.  the  proofs  of  the  strong  consistency  of  q  given 

2  2 

by  Zhao,  Krishnaiah  and  Bai  for  the  cases  of  unknown  o  and  a  =  1  are  valid  under  the 
following  weak  conditions: 

(i)  {x(t.,  i  =  1,2._}  is  a  stationary  and  ^  -  mixing  sample  sequence  with  ^  (3.11) 
decreasing  and  I  ^  (j)  < 

j-1 

•  244 

(ii)  E(x(t  I)  =  Cl  6({x'(t  lx(t  1}  )  <  •  f or  some  e  >  0  and  i  »  l.2,~  (3.12) 

•  1  •  •  1  •  1 

If  x(t  )_,x(t )  are  not  distributed  as  independent  complex  multivariate  normal,  L  and  L*  are 

-  1  -  n  k  k 

not'  likelihood  ratio  test  statistics  for  and  H*  but  the  procedures  described  in  this 
section  can  be  still  applied  if  the  conditions  (3.1 1)  and  (3.12)  are  satisfied 


4.  DETECTION  OF  NUNBB^  OF  SIGNALS  IN  PRESBSICE  OF  SPATIALLY  COLORS)  NOISE 


Consider  the  model 

xM  =  AsW  +  (4.1) 

where  X  is  known  or  unknowa  Also.  x(t)  :  pxl.  A.  s(t)  qx1  and  n(t)  are  as  defined  in 

Section  3  with  the  following  difference.  The  covariance  matrix  of  n(t|  is  Z  .  that  is,  noise 

•  1 

is  spatially  correlated  Here,  the  covariance  matrix  of  xlt)  is  related  to  Z^  as 

Z^  =  A^,A  +  XZ,  (4.2) 

Then,  we  are  interested  in  finding  out  the  number  of  signals  q  which  is  unknowa  This 
problem  is  equivalent  to  the  problem  of  finding  the  number  of  eigenvalues  of  Z^S^  'which 
are  different  from  its  smallest  eigenvalue. 

Rao  (19831  considered  the  following  problem  which  is  very  closely  related  to  the 
above  problem  Let  S^  and  S^  be  distributed  independetly  as  central  Wishart  matrices  with 
n  and  n^  degrees  of  freedom  respectively  and  E(S)  =  nZ  for  i  s  1.2.  Also,  let  Z  =  T 

12  i  i  i  2 

+  XZ^  where  X  is  known  or  unknown  and  T  is  of  rank  q(<p)  and  is  nonnegative  definita 
Rao  (19831  derived  the  likelihood  ratio  test  and  modified  likelihood  ratio  test  for  the  rank 
of  r  when  X  is  unknown  and  known  respectively.  The  above  procedures,  with  minor 
modifications,  can  be  applied  to  the  case  when  and  are  distributed  as  complex 
Wishart  matrices.  So.  the  results  of  Rao  can  be  applied  for  testing  the  hypothesis  that  the 
number  of  signals  is  equal  to  a  specified  value.  We  will  now  describe  certain  model 
selection  procedures  proposed  by  Zhao,  Krishnaiah  and  Bai  (1985b)  for  determination  of 
the  number  of  signals. 

When  the  data  is  collected  at  widely  spaced  time  points  t  ^it  is  realistic  to 

1  1^2 

assume  that  x(t  ) _ x|t  ),  the  observations  on  the  sensors,  are  distributed  independently.  In 

-  1  - 

"2  _ 

this  case,  we  can  estimate  Z.  with  S.(n.  where  S  =  7x|t)x(t)’.  Now.  let  S  /n  denote 

2  22  2  *•  -  j-  j  11 

1*1 

an  unbiased  biased  estimate  of  Z^  obtained  from  an  independent  set  of  n^  observations 

which  consist  of  noise  only  and  no  signals.  In  this  case  S  and  S.  are  distributed 

12 

independently  as  complex  Wishart  matrices  with  n^  and  n^  degrees  of  freedom  and  E(S.)  ^ 

nZ.  for  i  =  1,2  where  the  relationship  between  Z  and  Z.  is  given  by  (3.2)  with  Z 
II  12  p 


*.  %  %.  •% 


9 


instead  of  o\  Now.  let  X,  >  ...  >  X  denote  the  eigenvalues  of  Z^Z.'and  let  £  > 

1  p  *  2  1  1 

>  5  denote  the  eigenvalues  of  S.S~  n  /n  .  Also  let 

p  ^  2  112 


H;X>X  =.-=X=1. 

q  q  q+i  p 


(4.3) 


Zhao.  Krishnaiah  and  Bai  (1985bl  showed  that  the  LRT  statistic  for  is  given  by 


P 

q  !■  1 n  n  I  i 


(4.4) 


where  n  =  n^+  n^.  =  n^/n.  =  n^/n  and  t  denotes  the  number  of  fi's  greater  than 

one.  Now.  let 


EDC*(k.C  )  =  -logL*  +  Y*lc.p)C 

n  k  n 


(4.5) 


where  y  (k.p)  =  k(2p-k)  and  C  satistfies  the  following  conditions: 

n 


(i)  lim  1C  /n)  =  0 


(4.6) 


(ii)  lim  1C  /loglogn)  -  ». 
ru^  " 


(4.7) 


An  estimate  of  q  of  q  is  given  by 


roc  lq,C  )  =  min{roC  lO.C  lp-1,C  )}. 

n  n  n 


(4.8) 


Zhao.  Krishnaiah  and  Bai  (1985b)  proposed  the  above  estimate  and  proved  its  strong 
consistency. 


Next  let  H  :  X  >  X  , 

q  q  q+1 


~  X  =  X  where  X  is  unknown,  in  this  case,  the 

p 


likelihood  ratio  test  statistic  is  known  (see  Rao  (1983))  to  be 


L  =  '  E,  C((n,5+nX, 
q  l■q•*•1  2  I  I  kO 

A 

where  X^^  satisfies  the  equation 


(4.9) 


A  A 


Now,  lot 


=  I  {X  /(a  X,„  +  B  «.)>• 

i-iTM  ^’0  ' 


(4.10) 


EDC(k,C  )  =  -logL^  +  C 

n  ic  n 


(4.11) 


whore  Y^>P)  ~  l((2p-k)+1.  Then,  an  estimate  of  q  of  q  is  given  by 

EDC(q.C  )  =  min {HXiO.C  ).-.EDC(p- 1  ,C  )} .  (4.12) 

n  n  n 

Zhao,  Krishnaiah  and  Bai  (1985b)  proposed  the  above  estimate  and  established  its  strong 
consistency.  The  property  of  strong  consistency  holds  good  even  if  the  assumptions  of 
independence  and  normality  are  violated  provided  the  conditions  (4.6)  and  (4.7)  are  satisfied 
and  (3.11)  and  (3.12)  are  true  for  the  two  sets  of  observations  used  to  estimate  and 


We  will  now  discuss  the  method  proposed  by  Paulraj  and  Kailath  (1986)  for  the 
determination  of  the  number  of  signals  when  the  noise  is  spatially  colored.  They 
considered  the  situations  when  two  estimates  of  the  array  covariances  are  available  by 
displacing  the  array.  The  two-measurement  model  considered  by  them  is  given  by 


x(t)  =  As(t)  +  fKt)  i  =  1,2 


(4.13) 


where  x(t)  and  s.(t)  are  respectively  the  observation  vector  and  vector  of  wavefront  of  i-th 
measurement  The  covariance  matrices  of  the  observation  vectors  for  the  measurements  is 
given  by 


=  Ai»A  +  £  . 

2i  II  I  1 


(4.14) 


Also,  £.  is  estimated  with  S  /n  where 

2i  2i 


nS,.  =  I  x(t  )x(t  r. 

2i  .  •*  *1  j-i  j 
J*1 


In  addition,  let  P  =  £,.  -  £,^  and  R  =  (l/nMS^.  -  S).  Now,  let  X  >  _  >  X  denote 

21  22  21  22  1  p 

the  eigenvalues  of  P  and  let  >  ._  >  denote  the  eigenvalues  of  R  if  the  columns  of 
[A^  A^]  are  distinct  and  p^2q,  then 


X  =  -  =  X  =0. 
q+ 1  p-q 


If  i|»^  and  are  positive  definite  and  ^  ~  1(1,.  *>  nonsingular  and  A,  -  A  ,  then 


(4.15) 


TAing  advantage  of  (4.15)  and  (4.16),  Paulraj  and  Kailatti  proposed  using  nnethods  similar  to 

those  used  in  Wax  and  Kailath  (1984).  But  unfortunately,  the  distribution  of  -  S.  is 

22  21 

very  complicated  since  it  is  the  difference  between  two  complex  Wishart  matrices.  So, 
we  feel  that  the  derivation  of  the  likelihood  ratio  test  statistics  for  testing  the  hypothesis 
that  some  of  the  eigenvalues  are  zero  and  the  establishment  of  the  consistency  of  the 
information  theoretic  procedures  are  difficult  and  may  not  be  feasible.  When  »  A^, 

one  may  use  some  suitable  function  of  £  to  test  the  hypothesis  that  X  =  -  = 

k+l  p  k+l 

X^  =  0.  One  possible  test  statistic  is  But  unfortunately,  the  distributions  of  test 

statistics  based  upon  £  ,_,£  involve  X  ,~,X  as  nuisance  parameters.  But  asymptotic 

k+1p  Ik 

distributions  of  the  test  statistics  may  not  involve  nuisance  parameters.  Remarks  similar  to 
the  above  may  be  made  for  testing  the  hypothesis  that  X  =  ~  =  X  >0  whan  tk  and 

k+1  p-q  1 

4*2  o'a  positive  definite  and  the  columns  of  [A^  A^]  are  distinct 


5.  DETECTION  OF  NUMBER  OF  SIGNALS  HM  DECENTRALIZB}  PROCESSMG 

In  some  situations,  the  arrays  of  sensors  are  located  in  locations  which  are  wide 

apart  In  these  cases,  the  model  can  be  written  as 

X  W  »  A^s(t)  +  n  (t)  (5.1) 

for  h  3  1,2,..,M  where  x  (t)  :  p  x1  denotes  the  vector  of  observation  at  t-th  time  point 

-h  h 

using  h-th  group  of  arrays,  s(t)  :  qx  1  is  the  signal  vector,  n  (t)  is  the  noise  vector  on  h-th 
subarray  and  A  :  p  xq  is  a  matrix  whose  elements  depend  upon  the  parameters  of  the 

h  h 

signals.  We  also  assume  that  s(t)  and  n  (t)  are  distributed  independently  as  complex 

•  •11 

multivariate  normal  with  E(s(t))  «  Q  E(n  M)  =  (X  E(s(t)s(9k  -  and  E(n  (tlnlb')  =  ah.  We  also 

-  —  •n 

assume  that  x  (t)..-,x  (t)  are  distributed  independent  of  each  other.  The  covariance  matrix 
•  1  •M 

of  x^  is  given  by 


12 


Now  let  S  =  J  X  It  )x  It ).  AI$o,  let  £  >  -  >  £  ,  denote  the  eigenvalues  of  S  /a 

2h  J“1  h  j  h  j  1h  p  h 

h 

Under  the  model  15.1)  Wax  and  Kailath  11985)  proposed  to  estimate  q,  the  number  of 

sources,  using  AlC  and  MDL  criteria  The  same  method  can  be  used  even  if  the  noise 

variance  differs  from  subarray  to  subarray.  We  now  propose  a  method  of  estimation  of 

2 

the  number  of  signals  using  the  EDC  criterion  when  the  covariance  matrix  of  Mt)  is  o^. 

EDC(C^;p,l  =  -I09l.„  *  (5.3) 


where  C  satisfies  the  conditions 

nh 


li)  lim 


{C  /n}  =  0 


lii)  lim 
n.*.* 


{C  /loglogn}  =  «. 

nh 


(5.4) 


Also,  Ylk,p J  *  M2p -k+1)+  1,  and 

h  h 


In  addition,  let 


i-k+l 


M 

EDCIW  =  I  HX:(k,C  ;p  ) 

Nn  n 

h-1 

A 

Then,  an  estimate  q  of  q  is  given  by 


(5.5) 


(5.6) 


EOaq)  =  min{H)ClO),_,BXlp-1))  (5.7) 

where  p  =  minip  „-,p  ).  We  may  choose  C  s  to  be  equal  to  C  .  The  above  procedure 
may  be  applied  even  if  .  Alternatively,  we  can  use  the  procedure  described 

t 

below.  Let  denote  the  likelihood  ratio  test  statistics  for  testing  the  hypothesis  that 

X  *»  =  X  »o^h  =  1,2,_,Nl 
n 

In  this  case 


.'.V. 


■  .1 


.1 


L  =  {  n  n  a"/(  I  X  i  o 

k  h*!  i^k+l  ih  h»1  i*!^!  ih  0 


nip  -kMl 


where  p  =  p  +  ...  -t-  p  .  Now,  let 
0  1  M 


(5.8) 


(5,9) 


B)C(k,C  I  =  -logL^  +  C  y(k.p.) 

n  k  n  0 


where  vlk.p  )  =  M2p  -kM+M)  +  M  and  C  satisfies  the  condition 

'  0  0  n 


(i)  lim  <C  /n)  =  0 
n^«®  ” 

liil  lim  (C  /loglogni  =  ». 
n.^*  ” 

A  A 

Then,  the  number  of  signals  q  can  be  estimated  with  q  where  q  is  given  by 

HiClaC  )  =  min{H>ClO,C  ).-.,B3C(p-1,C  )}  (5.10) 

n  n  n 

where  p  =  minip  ,...,p  )  as  before.  Here,  we  note  that  x  (t),_,x  It)  are  not  distributed 
independent  of  each  other.  If  we  assume  that  n  |t)_,n  It)  are  distributed  independent  of 

each  other,  then  E{x|t)xW'l  =  AAaT’  So.  L  is  not  the  likelihood  ratio  tost  statistic  unless 

-i  -j  1  k 

A  i(»A.  =  0  for  i  ^j. 

Next  let  the  covariance  matrix  of  n  It)  be  and  otfier  assumptions  on  the  model 

•li  h 

15.11  remain  the  same  Also,  let  Z  denote  an  unbiased  estimate  of  Z,.  obtained  from  an 
independent  sample  of  observations  which  consist  of  noise  without  signal  We  assume  that 
the  size  of  this  sample  is  n,  and  the  observations  in  this  sample  are  independent  We  now 
propose  an  EDC  procedure  for  the  determination  of  the  number  of  signals  in  this  case. 

Lot  X  >  -  >  X  denote  the  eigenvalues  of  Z_Z,J  and  let  6  >  „  >  6 

Ih  D  h  2h  ih  Ih  p  h 


denote  the  eigenvalues  of  Z,^Z,J.  If  the  number  of  signals  is  q.  then 

2n  m 


=  ^  K  =  ^ 

V'’ 


(5.11) 


for  h  =  1.2.«,M.  Now.  let 


(5.12) 


where  y.  Ik,p  )  = 


w  M  ,~n/2_~nb  /2 

L  =  ^  Ka  +B  fi  j  5  " 

kh  l■1+mln<k,t  n  n  ih 
h 

A  A.  y 

and  T  denotes  the  number  of  eigenvalues  of  vvhich  are  greater 

h  2n  In 

C  satisfies  the  following  conditions: 

nh 


(i)  Urn  {C  /n)  =  0 

no.  ** 


(5.13) 


than  one  Also. 


(ii)  lim  {C  /ioglogn}  =  «. 

n  M  *** 


(5.14) 


Here,  we  note  that  L  denotes  the  likelihood  ratio  test  statistic  for  testing  the  hypothesis 


that  X  =  ~  =  X  =1-  Now. 
q+1.h  p  .h 

h 


EDCCk)  =  Z  EDC(k.C  ;p J. 

h«  t  nil  h 


Then,  an  estimate  q  of  q  is  given  by 


HXriql  =  miniaXiOUHJCip-l)} 


(5.15) 


(5.16) 


where  p  -  min(p^ . p^)  as  before.  By  following  the  same  lines  as  in  Zhao.  Bai  and 

Krishnaiah  (1985b).  the  consistency  of  the  above  procedure  can  be  established 


Now.  let  Z  s  ...  =  Z  and  let  L  denote  the  likelihood  ratio  test  statistic  for 

1 


testing  the  hypothesis 


In  addition,  let 


X  =~=X  =lforh  =  1.2_.M. 

k+1.ti  p.h 


EDCIk.C  )  =  -logL^  +  C  yfcl 

n  fc  n 


where  denotes  the  number  of  free  parameters  that  have  to  be  estimated  and 
satisfies  the  following  conditions; 


W  Km  {C  /n)  =  0 
a**  " 


(ii)  lim  {C  /Ioglogn)  =  ». 

A  A 

Then,  we  can  estimate  q  with  q  where  q  is  given  by 

aXldC  )  =  min{B3C|0.C  ).._.HX:ip-1,C  ) 

n  n  n 
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